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Abstract. Field-angle dependent specific heat measurement has been done on 
the heavy-fermion superconductor CeCoIns down to ~ 0.29 K, in a magnetic field 
rotating in the tetragonal c-plane. A clear fourfold angular oscillation is observed 
in the specific heat with the minima (maxima) occurring along the [100] ([110]) 
directions. Oscillation persists down to low fields H <C H c 2, thus directly proving 
the existence of gap nodes. The results indicate that the superconducting gap 
symmetry is most probably of d xy type. 
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Recently a new class of heavy-fermion superconductors CeTIn 5 (T = Rh, Ir, and 
Co) has been discovered. While CeRhIn 5 P is a pressure-induced superconductor, 
Celrln 5 [2] and CeCoIn 5 [3] show superconductivity at an ambient pressure at 0.4 
and 2.3 K, respectively. Among these, unique properties of CeCoIns have attracted 
much attentions in these years. Various experiments such as specific heat 0] 
thermal conductivity 0] and NMR T\ measurements |H] have revealed that CeCoIns 
is an unconventional superconductor with line nodes in the gap. Together with the 
suppression of spin susceptibility below T c Cj , this compound is identified as a d- 
wave superconductor. Recent rf penetration depth measurement |Hj and the flux line 
lattice imaging study by small-angle neutron scattering |§] also seem to be consistent 
with the existence of line nodes running along the c-axis. 

In unconventional superconductors, identification of the gap-node structure is of 
fundamental importance in understanding the pairing mechanism. Regarding this 
issue, it has recently been pointed out that the zero-energy density of states (ZEDOS) 
in the superconducting mixed state exhibits a characteristic oscillation with respect to 
the angle between H and the nodal direction |10l 1111 [T2] . An intuitive explanation of 
this effect employs a Doppler shift of the quasiparticle (QP) energy spectrum due to 
the local supercurrent flow [101 II II IT2"lll3| . More quantitative analysis of the effect has 
been done recently by Miranovic et al |14| within the quasiclassical formalism, which 
incorporates also the contribution of QPs in the core states of vortices. Experimentally, 
the angular dependent ZEDOS can be probed by the thermal conductivity [T3HTB1ITT) 
or the specific heat JB] measurements in rotating magnetic field at low T. As 
for CeCoIn5, the angle-resolved thermal conductivity (k(9)) measurement ^J] has 
revealed a clear fourfold oscillation in a magnetic field rotating within the c-plane. 
It was argued that the line nodes exist along the [110] directions (d 3 .2_ y 2 type gap 
symmetry). 

The interpretation of the thermal conductivity data is, however, necessarily 
involved because k(9) is proportional to the specific heat C{0) as well as to the QP 
scattering time r(9). The latter also depends on the field orientation but with opposite 
angular oscillation amplitude to the former. Which of these two predominantly 
contributes to n(0) is a subtle question, making it difficult to identify the gap node 
direction in some cases. It is therefore of vital importance to directly observe ZEDOS 
by the angle-resolved specific heat measurement. In this Letter, we have examined 
the C(H,9) of CeCoIns in a magnetic field rotating within the c-plane. We observed 
a clear fourfold oscillation in C (H, 9) , and argue that the oscillation originates from 
the nodal gap most probably of d xy symmetry. 

A single crystal of CeCoIns was grown by the so-called self- flux method |19| . 
The sample was cut into a thin plate (~ 2.3 x 2.0 x 0.7 mm 3 ) with the largest plane 
oriented perpendicular to the c-axis. Field-angle-dependent specific heat measurement 
was done by a standard adiabatic heat-pulse technique using a 3 He refrigerator. A 
transverse magnetic field was generated by a split-pair superconducting magnet and 
applied parallel to the c-plane of the sample mounted on a quartz platform. The 
3 He refrigerator was mounted on a mechanical rotating stage driven by a computer 
controlled stepping motor at the top of the Dewar, by which a quasi-continuous change 
of the field direction with the minimum step of 0.04° could be made. The angular 
dependent specific heat data were collected at an interval of 2° in the range ±105° 
with respect to the a-axis. At the temperature of 380 mK, the 2°-rotation of the 
refrigerator caused a heating of the sample by ~ 100 mK in a field of 5 T. In order 
to accurately evaluate C(H, 9) at constant temperature, each data point was taken 
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after the sample temperature relaxed to within 0.1% of the initial temperature. The 
details of the experimental set-up will be published elsewhere [201 • Field dependence 
of the specific heat up to 12.5 T was also measured by the relaxation method using a 
commercial calorimeter (PPMS, Quantum Design Co.). 

Figure shows the field dependence of the specific heat C(H) of CeCoIn 5 
obtained at T = 0.38 K for H || [100] and [110]. Throughout this paper, the 
nuclear-spin contribution C nuc has been subtracted from the data, assuming the 
form C* nuc = {A a + AxH^/T 2 with A a = 7.58 x 10~ 2 mJKmoh 1 and A 1 = 
6.90 x 10~ 2 mJ K mol -1 T~ 2 C(H) of CeCoIns in the superconducting mixed state 
is quite unusual. This behavior is very different from the one expected for ordinary s- 
wave superconductors in which C{H) linearly increases with H due to the contribution 
of QPs trapped in the vortex cores. In many of anisotropic superconductors, a power 
law dependence of C{H) oc H$ with (3 ~ 0.5 has been observed and is attributed to 
the nodal QP excitations. The negative curvature of C(H) in low fields seen in figure^ 
is thus consistent with the existence of gap nodes in CeCoIn5. However, the curvature 
of C(H) changes sign above 2 T and becomes strongly positive at higher fields, as 
shown in the inset. Similar behavior is previously reported for H || [001] [5], as well 
as in Sr 2 Ru04 for H || [110] Although the reason for this field dependence of 

the specific heat in CeCoIn 5 is unclear at the moment, it might be related to non- 
Fermi liquid behavior near the quantum critical point. With increasing H in the 
c-plane, we observed a small but discernible anisotropy in C(H) as shown in figure^ 
C(H || [100]) < C(H || [110]). 

The in-plane anisotropy of C (H) can be demonstrated more clearly by measuring 
its field-angular dependence, and some of the results are shown in figure First 
of all, we examined the contribution of the addenda (the lower trace of figure |^ a)), 
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Figure 1. Field dependence of the specific heat C(H) of CeCoIns for H \\ [100] 
(solid circles) and [110] (open circles), measured at T = 0.38 K. Inset: Overall 
field variation of C(H) up to 12.5 T for H \\ [100], measured at 0.5 K . 
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Figure 2. (a) Angular dependence of the total specific heat of CeCoIns (solid 
circles) mounted on the platform, in a field of 2 T rotated in the c-plane at 
T = 0.38 K. 9 = is the [100] direction. Open circles are the results without 
the sample, (b) Field-orientational dependence of the specific heat of CeCoIns 
for H = 2 T measured at T = 0.29 K and 1 K, and for H = 5 T measured at 
T = 0.38 K. Solid lines are the fit to C(H, 6)/T = (Co + Ch (1 + A 4 cos40))/T, 
where Co = 0.0361,0.0631 and 0.4042 Jmol^K- 1 for T = 0.29, 0.38 and 1 K, 
and A 4 = -0.0217 (2 T, 0.29 K), -0.0157 (5 T, 0.38 K) and -0.0061 (2 T, 1 K). 
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which is composed of the quartz platform, a thermometer (Ru02 chip resistor) and 
a heater. As can be seen in the figure, there is virtually no angular dependence 
in the addenda specific heat, implying that the field-orientational dependence of the 
thermometer is negligible. By contrast, there is a small but distinct fourfold oscillation 
in the raw specific heat data with the sample on the platform (the upper trace of 
figure EI a)). Comparing these two data, it is obvious that the observed angular 
oscillation is intrinsic to the sample. No appreciable twofold component is observed 
in the oscillation, indicating that the magnetic field is well oriented along the c-plane. 

The specific heat C(H, 9) /T of CeCoIn 5 is shown in figure^b) for H = 2 T (T = 
0.29 and 1 K) and 5 T (T = 0.38 K) as a function of the field angle 9 measured with 
respect to the a-axis. The fourfold oscillation is clearly seen at the base temperature of 
0.29 K but rapidly fades away at higher temperatures. C(H, 9) can be decomposed into 
a constant and field-angle-dependent terms: C(H, 9) = Co + C#(l + A{9)). Co is the 
zero- field term mainly due to thermally excited quasiparticles and phonons, whereas 
Gh and A(9) are field-dependent. The solid lines in figure |2Jb) are the fitting results 
assuming a simple form A{&) = A^ cos 49, by which we evaluated the amplitude of the 
fourfold angular oscillation. The sign of A4 is negative, with the minima of C(H, 9) 
occurring along the a-directions. In figure 01 we plotted the field dependence of the 
relative amplitude \A±\ as a function of H/H C 2- Most remarkably, |Aj decreases 
monotonously with H within the range of fields examined. The inset of figure [3] 
shows the temperature variation of the absolute amplitude \A±\ Ch/T measured at 
H = 2 T (H/H c2 = 0.175). The fourfold angular oscillation in C{H,6)/T rapidly 
diminishes with increasing T and vanishes above T c , implying that superconductivity 




Figure 3. Field dependence of the amplitude \A&\ of the fourfold oscillation at 
T = 0.38 K, plotted as a function of reduced field H/H c 2- The broken line is 
the calculated field variation for IA4I at T = assuming the d^^-type gap node. 
Inset: temperature dependence of the absolute amplitude \A4\Ch/T measured 
at H/H c 2 = 0.175. The solid line is a guide for the eyes. 
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is responsible for its occurrence. 

We now discuss the origin of the fourfold oscillation in C(H,9). First of all, 
we should pay attention to the in-plane anisotropy of H C 2, which is determined both 
by the gap topology and Fermi surface anisotropy; H C 2 || [100] is about 2.7% larger 
than H C 2 || [110] 19 . The in-plane anisotropy of H C 2 alone would also give rise to 
the angular dependence of ZEDOS. This is because the gap amplitude, which is a 
decreasing function of the reduced field H/H C 2, becomes field-angle dependent under 
fixed H . This effect certainly becomes important when H is near H C 2, and enhances 
C(H,8) for H || [110] along which the system is closest to the normal state; the sign 
of the fourfold amplitude is the same as what we observed. However, in the limit of 
low fields H <C H c2 the effect of the upper critical field can be disregarded. A crude 
estimate, assuming that C (H) in figure^can be scaled linearly with H C 2 (9) , shows that 
the ~ 3% anisotropy in H C 2 would give the fourfold term A4 of the order of 1% at 2 T. 
This definitely overestimates the effect of H C 2 anisotropy, but is still smaller than the 
observed one. Moreover, one expects the effect of the upper critical field to proliferate 
with increasing field. Instead, the oscillation amplitude is a decreasing function of H 
up to at least H ~ 0.5H C 2. The strong temperature variation of the fourfold amplitude 
of C(H)/T shown in the inset of the figure 01 is also hardly explained by the H C 2 in- 
plane anisotropy effect, because the upturn in C(H)/T does not change much in these 
temperature regions Since the effect of H C 2 anisotropy is incompatible with the 
data, we are led to conclude that the observed angular oscillation in C (H, 9) originates 
from the nodal structure, which has fourfold symmetry in the afo-plane. 

In the mixed state of anisotropic superconductors enhanced contribution to 
ZEDOS is coming from the QPs with momentum along the gap node direction. The 
enhancement depends on the relative position of the gap and magnetic field direction. 
Accordingly, ZEDOS becomes field-angle dependent, and exhibits a characteristic 
oscillation with respect to the angle between H and gap-node directions. This 
ZEDOS oscillation is directly reflected in the field-orientational dependence of the 
specific heat at low T, which takes minima (maxima) for H parallel to the nodal 
(antinodal) direction [TUJ H2 E3]. The clear 

minima of C(H,9) along [100], which 
persist down to low fields, indicate that there are gap nodes along these directions. 
Here we performed the microscopic calculation of the field-angle dependent ZEDOS 
of clean d-wave superconductor with spherical Fermi surface. We solved numerically 
and self-consistently the quasi-classical Eilenberger equations. This formalism is a 
good approximation as long as 3> 1 (kp being the Fermi wave number and £ the 
coherence length), the condition which is met in CeCoIns. It takes into account both 
the Doppler shift effect and the vortex core contribution on equal footing and without 
any further approximation. The result is presented with the broken line in figure 01 
Although the predicted curve gives slightly larger amplitude than the experimental 
one, it explains the overall field dependence of our experimental result reasonably 
well. The discrepancy between the predicted and observed amplitude may partly arise 
from the fact that the measurements have been done at finite T (T/T c ~ 0.17), as 
inferred from the strong temperature variation of the amplitude predicted in Ref. ^21 • 
The actual amplitude however gradually levels off below 0.5 K. Impurities, which are 
always present, may reduce the oscillation amplitude at lower T. The present data 
thus strongly indicate that the gap symmetry of CeCoIn 5 is most probably of d xy 
type. 

The present results agree with the k(9) measurement on CeCoIns ^7] in the point 
that the nodal structure has fourfold symmetry in the c-plane, but disagrees on the 



Letter to the Editor 



7 



location of the nodes. k(9) shows a fourfold oscillation with the minima along the 
[010] and [100] directions, which is the same oscillation behavior as our C(9) data. 
The authors of Ref. J7| assumed that this angular dependence of k(9) is dominated 
by the QP scattering time, which becomes largest (smallest) when the magnetic field 
is along the nodal (antinodal) directions, just opposite to the ZEDOS contribution. 
It is considered that the QP scattering effect is important at high temperatures, but 
the ZEDOS contribution becomes predominant with decreasing T . If the angular 
oscillation of k(6) had come from the QP scattering term, then its amplitude should 
change sign on cooling. The observed amplitude of the fourfold term in k(9) however 
continues to increase down to the lowest T of 0.35 K ^7]. Whether the amplitude 
changes sign at still lower temperatures or not would be an interesting issue, but 
our C(9) data implies that the ZEDOS effect is already predominant in n(9) at 
temperatures above 0.35 K. 

Present data also does not contradict the neutron scattering experiment by 
Eskildsen et al UJ, who find the square vortex lattice oriented along the [110] 
direction. The result of neutron scattering experiment was argued to support the 
picture of CeCoIns as being the d x 2_ y 2 superconductor. However, the structure and 
the orientation of the vortex lattice is determined by the combined effect of the Fermi 
surface anisotropy and the gap structure. Which effect will prevail is a subtle question 
and the answer depends on many factors: field value, temperature, degree of Fermi 
surface anisotropy etc. Thus, the orientation of the square vortex lattice in tetragonal 
crystals may not serve as a conclusive test for the positions of gap nodes. Model 
calculation by Nakai et al t 22] shows a variety of stable vortex lattice structures in 
tetragonal d xy -wa.ve superconductors. It is also shown that the low field square vortex 
lattice is oriented along [110] direction which is consistent with the neutron scattering 
experiment. It is interesting to point out that the square vortex lattice orientation 
should rotate by 7r/4 in high fields if the d xy is indeed correct identification of the gap 
function. 

In summary, we have performed the angle-resolved specific heat measurements on 
the heavy-fermion superconductor CeCoIns in a magnetic field rotating in the basal 
a6-plane. We observed a clear fourfold symmetry in C(H, 9) with the minima oriented 
along [100] directions, which comes from the field-angular oscillation of the zero-energy 
density of nodal quasiparticles. The results imply that the superconducting gap node 
of CeCoIn 5 is located most likely along the [100] and [010] directions, suggesting the 
symmetry to be of d xy type. 
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